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Abstract 

In this paper, which completes our earher short pubhcation [Phys. Rev. Lett. 84, 
511 (2000)], we study dynamics of a hard-core tracer particle (TP) performing a biased 
random walk in an adsorbed monolayer, composed of mobile hard-core particles under- 
going continuous exchanges with a vapor phase. In terms of an approximate approach, 
based on the decoupling of the third-order correlation functions, we obtain the density 
profiles of the monolayer particles around the TP and derive the force- velocity relation, 
determining the TP terminal velocity, Vtr, as the function of the magnitude of external 
bias and other system's parameters. Asymptotic forms of the monolayer particles density 
profiles at large separations from the TP, and behavior of Vtr in the limit of small external 
bias are found explicitly. 



PACS numbers: 05.40-a, 66.30.Lw, 68.45.Da 



1 Introduction. 



When a gas or vapor is brought into contact with a clean sohd surface, some portion of the gas 
particles becomes reversibly attached to the surface forming an adsorbed layer. Following the 
seminal work of Langmuir (see, e.g., Ref.|j|]), equilibrium properties of such adsorbed layers 
have been extensively studied and a significant number of important advancements have 
been made. In particular, subsequent analysis included more realistic forms of intermolecular 
interactions or allowed for a possibility of multilayer formation. As a result, different phase 
transformations have been predicted and different forms of adsorption isotherms have been 
deduced, which are well corroborated by available experimental data (see, e.g. Refs.|l]-^). 

Some work has been focused on understanding of random motion of individual molecules 
in adsorbed layers, which constitutes an important factor limiting their global dynamical 
behavior. For instance, transport processes control spreading rates of molecular films on solid 



surfaces ^, spontaneous or forced dewetting of monolayers [^-^ or island formation |10|. 
Here, some approximate analytical results have been obtained for both dynamics of an isolated 
adatom on a corrugated surface and collective diffusion, describing spreading of macroscopic 
density fluctuations in interacting adsorbates being in contact with the vapor phase 
p^ . Analysis of tracer diffusion in adsorbed layers, which provides a useful information 
about intrinsic friction of the adsorbates, pertains mostly to two-dimensional models with 
forbidden particles exchanges with the vapor (see, e.g. Refs.||-|2|. Tr acer diffusion in 
adsorbed monolayers undergoing continuous adsorption and desorption has been essentially 
less studied. A few available detailed studies are either devoted to an analysis of a somewhat 
artificial one-dimensional situation |21|, or, for two-dimensional adsorbates in contact with 
the vapor, focus on the limit when the tracer particle perform a totally directed random walk 



j22|. The results for the general situation with arbitrary bias has been only briefly presented 
in Ref.^l. 

In this paper, which completes our earlier short publication |23|, we present a detailed 
analysis of the properties of a biased tracer diffusion in two-dimensional adsorbed monolayers 
undergoing continuous exchanges with the vapor. The system we consider consists of (a) a 
solid substrate, which is modeled in a usual fashion as a regular square lattice of adsorbtion 
sites; (b) a monolayer of adsorbed, mobile hard-core particles in contact with a vapor and 
(c) a single hard-core tracer particle (TP). The monolayer particles move randomly along the 
lattice by performing symmetric hopping motion between the neighboring lattice sites, which 
process is constrained by mutual hard-core interactions, and may desorb from and adsorb 
onto the lattice from the vapor with some prescribed rates dependent on the vapor pressure, 
temperature and the interactions with the solid substrate. In contrast, the tracer particle is 
constrained to move on the two-dimensional lattice only, (i.e. it can not desorb to the vapor), 
and is subject to a constant external force E. Hence, the TP performs a biased random walk, 
constrained by the hard-core interactions with the monolayer particles, and always remains 
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within the monolayer, probing its frictional properties. 



In terms of an approximate approach of Ref.|24|, based on the decouphng of the third- 
order tracer-particle-particle correlation functions into the product of corresponding pairwise 
correlations, we determine the density profiles of the monolayer particles, as seen from the 
stationary moving tracer, and calculate analytically the terminal velocity Vtr attained by the 
tracer particle in the limit t ^ oo. We show that the monolayer particles distribution is 
strongly inhomogeneous: the local density of the monolayer particles in front of the tracer 
is higher than the average, which means that the monolayer particles tend to accumulate in 
front of the driven tracer, creating a sort of a "traffic jam", which impedes its motion. The 
condensed, "traffic jam" -like region vanishes as an exponential function of the distance from 
the tracer. The characteristic length and the amplitude of the density relaxation function are 
calculated explicitly. On the other hand, past the tracer the local density is lower than the 
average. Here, we find that depending on whether the number of particles in the monolayer 
is explicitly conserved or not, the local density past the tracer may tend to the average value 
at large separations from the tracer in a completely different fashion: In the non-conserved 
case the decay of the density is described by an exponential function, while for the conserved 
particles number case it shows an algebraic dependence on the distance, revealing in the latter 
case especially strong memory effects and strong correlations between the particle distribution 
in the monolayer and the tracer position. Further on, we find that the terminal velocity of 
the tracer particle depends explicitly on both the excess density in the "jammed" region in 
front of the tracer, as well as on the density in the depleted region past the tracer. We 
realize that both densities are themselves dependent on the magnitude of the tracer velocity, 
applied external force, as well as on the rate of the adsorption/desorption processes and 
on the rate at which the particles can diffuse away of the tracer, which results in effective 
non-linear coupling between Vtr and E. In consequence, in the general case (for arbitrary 
adsorption/desorption rates and arbitrary external force), Vtr can be found only implicitly, as 
the solution of a transcendental equation relating Vtr to the system parameters. This equation 
simplifies considerably in the limit of a vanishingly small external bias; in this case we obtain 
a linear force-velocity relation, akin to the so-called Stokes formula, which signifies that in 
this limit the frictional force exerted on the tracer particle by the host medium (the adsorbed 
monolayer) is viscous. The corresponding friction coefficient is also determined explicitly. 

We finally remark, that a qualitatively similar physical effect has been predicted recently 
for a different model system involving a charged particle moving at a constant speed a small 
distance above the surface of an incompressible, infinitely deep liquid. It has been shown 



in Refs.|25, 26|, that the interactions between the moving particle and the fiuid molecules 
induce an effective frictional force exerted on the particle, producing a local distortion of the 
liquid interface, - a bump, which travels together with the particle and increases effectively 
its mass. The mass of the bump, which is analogous to the jammed region appearing in our 
model, depends itself on the particle's velocity resulting in a non-linear coupling between the 
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medium-induced frictional force exerted on the particle and its velocity p5|, p^. 

The paper is structured as follows: In Section 2 we formulate the model and introduce 
basic notations. In Section 3 we write down the dynamical equations which govern the time 
evolution of the monolayer particles and of the tracer, and outline the decoupling proce- 
dure. Section 4 is devoted to the analytical solution of the decoupled discrete-space evolution 
equations in the limit i ^ oo in terms of the generating function approach. In Section 5 
we analyse the asymptotical behavior of the density profiles at large separations in front of 
and past the tracer particle. Next, in Section 6, we derive formal expressions describing the 
shape of the density profiles in the adsorbed monolayer and a force-velocity relation in the 
general case. Section 7 is devoted to the analysis of the general force- velocity relation in the 
limit of a vanishingly small external bias. Here, we derive an analog of the Stokes formula 
for biased diffusion in adsorbed monolayers undergoing continuous exchanges with the vapor 
phase, determine the friction coefficient and estimate, through the Einstein relation, the tracer 
diffusion coefficient. Finally, we conclude in Section 8 with a brief summary and discussion 
of our results. 



2 The model 

Consider a two-dimensional square lattice of adsorption sites of spacing a, which is brought in 
contact with a reservoir containing identic, electrically neutral particles (vapor phase) (Fig.l), 
maintained at a constant pressure. We suppose that (a) the particles may leave the reservoir 
and adsorb onto any vacant lattice site at a fixed rate f/r*; (b) the adsorbed particles may 
move randomly along the lattice by hopping at a rate l/4r* to any of four neighboring sites, 
which process is constrained by hard-core exclusion preventing multiple occupancy, and (c) 
the adsorbed particles may desorb from the lattice back to the reservoir at rate g/r*. Both 
/ and g are site and environment independent. 

To describe the time-dependent occupancy of lattice sites, we introduce the variable ?7(R), 
which may assume two values: 




1, if the site R is occupied by an adsorbed particle, 
0, if the site R is empty. 



Evidently, local 7?(R) can change its value due to adsorption, desorption and constrained 
random hopping events, and the total number of particles in the adsorbed monolayer is not 
explicitly conserved due to adsorption/desorption processes. On the other hand, the mean 
density ps of the adsorbate, ps =< ??(R) >, approaches as t ^ oo a constant value 



which relation represents the customary Langmuir adsorption isotherm . We note also that 
in the analysis of the stationary-state behavior, we can always tmm to the conserved particles 
number limit by setting / and g equal to zero and keeping their ratio fixed, i.e. supposing 
that f/g = Ps/{^ ~ Ps)- This limit will correspond to the model of biased tracer diffusion in 
a two-dimensional hard-core lattice gas with fixed particles density ps, and will allow us to 
check our analytical predictions against some already known results p5|-p0|]. 

Further on, at t = we introduce at the lattice origin an extra hard-core particle, whose 
motion we would like to follow here and whose position at time t we denote as Titr- This 
tracer particle - the TP, can be thought of as an external probe designed to measure the 
resistance offered by the monolayer particles to the external perturbance, or, in other words, 
to measure the intrinsic frictional properties of the adsorbate. 

Now, we stipulate that the TP is different from the adsorbed particles in two aspects: 
first, it can not desorb from the lattice and second, it experiences an action of some external 
driving force, which favors its jumps into a preferential direction. Physically, such a situation 
may be realized, for instance, if this only particle is charged and the system is subject to a 
uniform electric field E. We suppose here, for simplicity of exposition, that the external force 
E is oriented in the positive X-direction, i.e. E = {E, 0). 

More precisely, we define the TP dynamics as follows: We suppose that the TP, which 
occupies the site Utr at time t, waits an exponentially distributed time with meanQ r, and 
then attempts to hop onto one of four neighboring sites, Utr + Gj^, where are four unit 
vectors of the square lattice. In what follows we adopt the notation = {±1, ±2}, where itei 
(respectively, =be2) will denote itX (respectively, ±y) directions. Next, the jump direction 
is chosen according to the probablity p^, which is defined in a usual fashion as 



exp 


f(E-e,)" 




exp 


f(E-e^)^ 



where /? is the reciprocal temperature, (E • e) stands for the scalar product, the charge of 
the TP is set equal to unity and the sum with the subscript p, denotes summation over all 
possible orientations of the vector e^, that is ;U = {±1, ±2}. 

After the jump direction is chosen, the TP attempts to hop onto the target site. The hop 
is instantaneously fulfilled if the target site is vacant at this moment of time; otherwise, i.e., if 
the target site is occupied by any adsorbed particle, the jump is rejected and the TP remains 
at its position. 

^We suppose that in the general case this mean time t is different from the corresponding time r* associated 
with the monolayer particles dynamics. As a matter of fact, this should be the case merely because the TP- 
substrate interactions may be different from the particle-substrate ones. Varying r we can mimic different 
possible situations; in particular, r = corresponds to the case when the tracer simply slides along the 
substrate regardless of the surface corrugation. 
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3 Evolution equations 



Let r] = {ri(R)} denote the entire set of the occupation variables, which defines the instan- 
taneous configuration of the adsorbed particles at the lattice at time moment t. Next, let 
P(Rtr,^;i) denote the joint probability of finding at time t the TP at the site Rtr and all 
adsorbed particles in the configuration rj. Then, denoting as rj''''^ a configuration obtained 
from r] by the Kawasaki-type exchange of the occupation variables of two neighboring sites 
r and r -|- ejy, and as 77"" - a configuration obtained from the original r] by the replacement 
ri{r) ^ 1 — rj{r), which corresponds to the Glauber-type flip of the occupation variable due 
to the adsorption/desorption events, we have that the time evolution of the configuration 
probability P(Rtr,^/;t) obeys the following master equation: 

M r^Rtr-e^,Rtr 

+ ^YPf^{i^- vCRtr)) P(Rtr - e^, 77; t) - (1 - r/(Rtr + e^)) P(Rtr, t)} 
+ ^11 {(l-7?(r))P(Rtr,r/'-;t)-r?(r)P(Rtr,r/;t)} 

+ E {r?(r)P(Rtr,r?'-;t)-(l-r/(r))P(Rtr,7?;t)}, (3) 

where the subscript r 7^ Rtr under the summation symbol signifies that summation extends 
over all lattice sites except for the site occupied at this time moment by the TP. 

Now, the instantaneous velocity Vtr{t) of the TP can be obtained by multiplying both sides 
of Eq.@ by (Rtr • ei) and summing over all possible configurations (Rtr,??)- This results in 
the following exact equation determining the TP velocity: 

= l J](Rtr-ei) = ^{pi(l-A:(ei;t)) -p_i(l-A;(e_i;t))}, (4) 

where 

k{X; t)=Y^ r/(Rtr + A)P(Rtr, V, t) (5) 

is the probability of having at time t an adsorbed particle at position A, defined in the frame 
of reference moving with the TP. Evidently, k{X) can also be interpreted as the density profile 
in the adsorbed monolayer as seen from the moving TP. 

Equation (^) signifies that the instantaneous velocity of the TP is dependent on the 
monolayer particles density in the immediate vicinity of the tracer. If the monolayer is 
perfectly stirred, or, in other words, if k{X) = ps everywhere, (which implies immediate 
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decoupling of Rtr and rj) , one would obtain from Eq. a trivial mean- field result 

Vl^'> = {p^-p^,){l-ps)^, (6) 

which states that the only effect of the medium on the TP dynamics is that its jump time r 
gets merely renormalized by a factor {I — ps)~^ , which represents the inverse concentration of 
voids in the monolayer; note that (1 — Ps)/t defines simply the mean frequency of successful 
jump events. However, the situation appears to be more complicated and, as we proceed 
to show, k{\) is different from the equilibrium value ps everywhere, except for |A| — > oo. 
This means that the TP strongly perturbs the particles distribution in the monolayer - it is 
no longer uniform and some non-trivial stationary density profiles emerge. Moreover, k{\) 
appears to be dependent on the TP velocity which results ultimately in a non-linear coupling 
between Vtrif) and density profiles. 

Now, in order to calculate the instantaneous mean velocity of the TP we have to determine 
the mean particles density at the neighboring to the TP sites Rtr + ej-i, which requires, in 
turn, computation of the density profile k{X) for arbitrary A. The latter can be found from the 
master equation (^) by multiplying both sides by r/(Rtr) and performing the summation over 
all configurations (Rtr,??)- In doing so, we find the following set of equations (see Appendix 
A for more details): 

AT*dtk{\;t) = ^(V^-5A,e^V_^)fc(A;t)-4(/ + 5)A;(A;t) + 4/ 

+ ^ J] J]p/.(l-??(R.tr + e^))V^7?(Rtr + A)P(Rtr,r/;t), (7) 

Rtr,»? M 

where denotes the ascending finite difference operator of the form 

V^/(A)^/(A + e^)-/(A), (8) 

and 

J 1, if the site r = r', 

<^>r,r' = \ 

I 0, otherwise. 

The Kroneker-delta term (^A,e^ signifies that the evolution of the pair correlations, Eq.(|^, 
proceeds differently at large separations and at the immediate vicinity of the TP, because 
of its asymmetric hopping rules, Eq.(^) (see for more details the points (a) and (b) in the 
Appendix A) . 

Note next that the contribution in the second line in Eq. ^ , associated with the TP biased 
diffusion, is non-linear with respect to the occupation numbers such that the pair correlation 
function gets effectively coupled to the evolution of the third-order correlations of the form 

T(A, e^; t) = ^(I^tr + A)7?(Rtr + e^)P(Rtr, r?; t). (9) 
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That is, Eq.(|^ is not closed with respect to the pair correlations but rather represents a 
first equation in the infinite hierachy of coupled equations for higher-order correlation func- 
tions. One faces, therefore, the problem of solving an infinite hierarchy of coupled differential 
equations and needs to resort to an approximate closure scheme. 

We proceed along the lines suggested in Ref . 1 23 ] and apply the simplest non-trivial closure 



approximation, based on the decoupling of the third-order correlation functions into the 
product of pair correlations. More precisely, we assume that for A 7^ the third-order 
correlation can be written down in the following form 



r/(Rtr + A)r/(Rtr + e^)P(Rtr, r?; t) 



^ r?(Rtr + A)P(Rtr,r/;t) r?(Rtr + e^)P(Rtr, r?; t) 

\Rtr,TI j \R.tr,»7 

= A;(A;t)A;(e^;t), (10) 

i.e. the joint probability of having at time moment t one adsorbed particle in the immediate 
vicinity of the TP, at position e^, and another particle at position A, is represented as the 
product of the corresponding pairwise probabilities. We hasten to remark that the approx- 
imate closure of the evolution equations in Eq.([lO|) has been already employed for studying 
related models of biased tracer diffusion in hard-core lattices gases and has been shown to 
provide quite an accurate description of both the dynamical and stationary-state behavior. 
The decoupling in Eq.(|lO|) has been first introduced in Ref.p^ to determine the properties of 
a driven tracer diffusion in a one-dimensional hard-core lattice gas with a conserved number 
of particles, i.e. without an exchange of particles with the reservoir. Extensive numerical 
simulations performed in Ref.|24| have demonstrated that such a decoupling provides quite 



a plausible approximation for the model under study. Moreover, rigorous probabilistic anal- 
ysis of Ref. [p7| has shown that for this model the results based on the decoupling scheme 
in Eq.([TO|) are exact. Furthermore, the same closure procedure has been recently applied to 
describe spreading kinetics of a hard-core lattice gas from a reservoir attached to one of the 
lattice sites B. Again, a very good agreement between the analytical predictions and numeri- 



cal results has been found. Next, the decoupling in Eq.(lO) has been used in a recent analysis 
of a biased tracer dynamics in a one-dimensional model of adsorbed monolayer in contact 
with a vapor phase p^, i.e. a one-dimensional version of the model to be studied here. Also 



in this case an excellent agreement has been observed between the analytical predictions and 



Monte Carlo simulations data [21|. Besides, as we have already mentioned in Ref. p3| , the 



closure of the hierarchy of the evolution equations in Eq.(|T^) allows us to reproduce in the 
limit f,g = and f /g = const (conserved particles number limit) the results of Refs.|16| and 



|17|, which are known (see e.g. Ref. p5| ) to provide a very good approximation for the tracer 



diffusion coefficient in two-dimensional hard-core lattice gases with arbitrary particle density. 
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Using the approximation in Eq.(|To|), we can rewrite Eq.(^) in the fohowing closed form 

4T*dtk{X;t) = Lk{X;t) +4f, (11) 

which holds for all A, except for A = {0, itei, ±62}. One the other hand, for these special 
sites A = e^, with 1/ = {±1, ±2} we find 

4T*dtkie^; t) = {L + A^)k{e^; t) + 4/, (12) 

where L is the operator 

L^^A^V^-4(/ + r7), (13) 

and the coefficients are defined by 

A^{t) = 1 + ^P^(l - k{e^; t)). (14) 
Now, several comments about Eqs.(p!l|) and (Il2|) are in order. First of all, let us note 



that Eq.(12) represents, from the mathematical point of view, the boundary conditions for 
the general evolution equation (lO) , imposed on the sites in the immediate vicinity of the TP. 



As we have noticed already, Eqs.(p^) and (12) have a different functional form since in the 
immediate vicinity of the TP its asymmetric hopping rules perturb essentially the monolayer 
particles dynamics. 

Next, Eqs.(|TT|) and (|l^ ) possess some intrinsic symmetries and hence the number of in- 
dependent parameters can be reduced. Namely, reversing the field, i.e. changing E — > — E, 
leads to the mere replacement of k{ei; t) by /c(e_i; t) but does not affect k{eu] t) with 1/ = ±2, 
which implies that 

A;(ei;t)(-E) = A;(e_i;t)(E), and fc(e^; t)(-E) = A;(e^; t)(E) for v = ±2, (15) 
Besides, since the transition probabilities in Eq.(^) obey p2 = one evidently has that 

/c(e2;t) = /c(e_2;t), (16) 

and hence, by symmetry, 

A2{t) = A^2{t) (17) 

which somewhat simplifies Eqs.(|Tl|) and (|l^). 

Lastly, we note that despite the fact that using the decoupling scheme in Eq.(|lO|) we 
effectively close the system of equations on the level of the pair correlations, solution of 
Eqs.dnJ) and (|l2|) still poses serious technical difficulties: Namely, these equations are non- 
linear with respect to the TP velocity, which enters the gradient term on the rhs of the 
evolution equations for the pair correlation, and does depend itself on the values of the 
monolayer particles densities in the immediate vicinity of the TP. Below we discuss a solution 
to this non-linear problem, focusing on the limit t 00. 
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4 Stationary solution of the evolution equations 



We turn to the limit t ^ oo and suppose that both the density profiles and stationary velocity 
of the TP have non-trivial stationary values 

k{X) = lim k{X;t), Vtr = lim Vtr{t), and An = lim AJt) (18) 

Define next the local deviations of A;(A) from the unperturbed density as 

h{X) = k{X) - ps (19) 

Choosing that h(0) = 0, we obtain then the following fundamental system of equations: 

Lh{X) = 0, (20) 

which holds for all A except for A = {0, e-ti, e-i-2}, while for the special sites adjacent to the 
TP, i.e. for A = {0, e-ti, e-1-2}) one has 

{L + A^)h{e^) + p,{A^ - A^^) = 0, (21) 

Equations ( pO[ ) and ( ^T[ ) determine the deviation from the unperturbed density ps in the 
stationary state. Note also that in virtue of the symmetry relations in Eqs.(p^) and (|l7|), 
h{e2) = h{e^2) and A2 = A^2- 

Now, our general approach to solution of coupled non-linear Eqs.(^),(^) and ( |2l| ) is as 
follows: We first derive a general solution of these equations supposing that Vtr is a given 
parameter, or, in other words, assuming that the coefficients A^, entering Eqs.(pO|) and ( ^l]) 
are known. In doing so, we obtain h{X) in the parametrized form 

h{X) = hiX;A±i,A2). (22) 



Then, substituting into Eq.(22) particular values A = {e-ti, e-1-2} and making use of the 
definition of A^ in Eq.(|l^, we find a system of three linear equations with three unknown 
coefficients of the form 

Au = l + ^Vu{l -Ps- /i(e^; A±i,^2)), (23) 

where v = {±1,2}, which will allow us to define all A^, (and hence, all h{e^)). Finally, 
substituting the results into Eq.(^), which can be written down in terms of Ay as 

^t'- = 4^(^1-^-1)' (24) 
we will arrive at a closed-form equation determining implicitly the stationary velocity. 
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4.1 Generating function for the stationary particle density profiles in the 
monolayer. 



Equations (]20|) and (21) can be most conveniently solved by introducing the generating func- 



tion for the density profiles of the form 

H{WI,W2)^ E Kun,w^i'w^^\ (25) 



-OO +00 



ni=— oo n2=— OO 



where /ini,n2 = h{nie\ + 71-262). Multiplying both sides of Eqs.(pO[) and (21) by w^^w'^'^ , and 
performing summations over ni and 722 we find that H{wi, W2) is given explicitly by 

H{wi,W2) = - K{wi,W2) ^Aiw^^ + A^iwi + A2{w2 + W2^) - , (26) 

where 

a^^^, + 4(/ + 5), (27) 

and 

K{wx,W2) = E^-Kh " ^^^^""^^ + P'^^^ ~ - ^r')- (28) 



Equations ( pq ) to ( p8D determine the generation function for the density profiles exactly, and 
the latter can be explicitly obtained via standard inversion formulae. 

4.2 Integral characteristic of the density profiles 

As we have already remarked, the presence of the driven TP induces an inhomogeneous density 
distribution in the monolayer. One can thus pose a natural question whether equilibrium 
between adsorption and desorption processes gets shifted due to such a perturbancy, i.e. 
whether the equilibrium density in the monolayer is different from that given by Eq.(|). The 
answer is trivially "no" in the case when the particles number is explicitly conserved, but in 
the general case with arbitrary / and g this is not at all evident: similarly to the behavior in 



one-dimensional system |21|, one expects that also in two-dimensions the density profiles are 
asymmetric as seen from the stationary moving TP and are characterized by a condensed, 
"traffic-jam" -like region in front of and a depleted region past the TP. One anticipates then 
that the desorption events are favored in front of the TP, while the adsorption events are 
evidently suppressed by the excess density. On the other hand, past the TP desorption is 
diminished due to the particles depletion while adsorption may proceed more readily due to 
the same reason. It is thus not at all clear a 'priori whether these two effects can compensate 
each other exactly, in view of a possible asymmetry of the density profiles, as it happens in 



the one-dimensional model |21|. 
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For this purpose, we study the behavior of the integral deviation of the density from 
the equihbrium value ps, i.e. 

+00 +00 

^= (29) 

ni=— 00 n2=— 00 

which can be computed straightforwardly from Eqs.(p6|) and ( pSj ) by setting both wi and W2 
equal to unity. Noticing that = l,tt;2 = 1) = 0, and that 2^2—0 = — 4(/+(7), 

i.e. is strictly negative as soon as adsorption/ desorption processes are present, we obtain 
then that Vt is stricly equal to 0. This implies, in turn, that the perturbancy of the density 
distribution in the monolayer created by the driven TP does not shift the global balance 
between the adsorption and desorption events. An analogous result has been obtained for the 
one-dimensional problem in Ref. ||2l|| . 



4.3 Inversion of the generating function with respect to the "symmetric" 
coordinate W2- 

Inversion of the generating function with respect to W2 can be readily performed by expanding 
H{wi,W2) into the series in powers of W2- To do this, we first rewrite H{wi,W2) as 

+00 ^ 
H{wi,W2) = a"^/C(i(Ji,'u;2)^a~*(Aiu;]"^ +yl2('W2 + 'w^^)) . (30) 

i=0 

Then, the sum on the right-hand-side of Eq.(^0|) can be further expanded as 

+00 



i=0 



w 



2 ' 



n2=— 00 j=0 j=0 



(31) 

where ("") stands for the binomial coefficients. Next, gathering the terms with the same power 



of W2 on the two sides of Eq. ( pOD and using Eq. (31), we have that 

+00 

Y V,n2<' = a"^{-X]^'^^(*^'^) + (^i^(^i) + '°^(^i-^-i^J^i 

ni=— 00 V 

+ (A_i/i(e_i) - ps{Ai - A^i))w^^]F^^\wi,\n2\) 

+ a-^A2h{B2)(F^'^\wu\n2-l\)+F^'^\wi,\n2 + l\)), (32) 

where 

F(^) («.!, n2) ^ g g + n2 + 2 A fn2 + 2 A ^n2+2, 1 + _ (33) 
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As the next step, we resum the series entering F^'^\wi,n2)- Substituting to Eq.(^ 
expUcit expressions for the binomial coefficients and for the Gamma function, we have 

n2 +00 , / /I \ 2i -, 



^ (^) E{(^) 



j=0 

00 

X / e-*t"2+2i 




a J r(i + i)r(n2+i + i) 



which yields 



X , 



Next, using the identity 



i=o V J / 

we obtain 

'^=\ , (37) 

where y = ^2/(0 — (Aiz/;^"^ + A^iwi)). Eqs.(|32[) to (^) define the inverted with respect to 
W2 generating function. 

4.4 Inversion of the generating function with respect to the "asymmetric" 
coordinate wi. 

Inversion of H{wi,W2) with respect to wi can be performed along exactly the same lines as 
we did it in the previous subsection for the variable W2- That is, by expanding F'^'^\wi^n2) 
in Eq. ( p7| ) into a series in powers of wi : 

+00 +00 +00 

ni=-oo fc=0 j=0 

\ni\ +n2 + 2j + 2k\ {n2 + 2j\ /|ni| + 2k'^ 



n2 + 2j J \ j J \ k 



A'^^+^^{A^A_i)''w'l'y (38) 



and then identifying the coefficients in this expansion with the analogous coefficients in the 
expansion in powers of the variable wi of 

hnun2 = J^^^/i(e^)V_^F„,,„2 - psiAi - A_i)(Vi - V-OK^.n^}. (39) 
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with 



fc=o i=o 



Then, using an integral identity 



nil + \n2\ + 2j + 2k\ /|n2| + 2j\ fni + 2k 
\n2\ + 2j J\ j J\ k 

^-t^\ni\ + \n2\+2k+2j+l 



T{k + i)r(|ni| + k + i)r(j + i)r(|n2| + i + 1) 

(41) 



we cast F„i,n2 into the form 

ni/2 



(42) 



where I„(2;) stands for the modified Bessel, defined as 



I„(z) = — /" de cos{n9)exp{zcos{9)) (43) 

271" J-TT 

It is worth-while to mention now that i^ni,n2 has an interesting physical interpretation. To 
illustrate it, we rewrite the integral involved in Eq.(^) using the definition in Eq.(|4^) as 

f-oo rn rn ^-{kini+k2n2) 

e-%,[2a-^y^A^itjln2i2a-^A2t)dt = j^ J J i _ ^x{k) ^^^^^2, (44) 

where 

„ _i / r-. — -. . , 1 , /, X COs(/ci) + A2 C0S(/C2) 

i = 2a i(V^i^_i+A2), and A(k) = nA^ , a ' (^^) 

One recognizes then (cf . ref . ||2^ ) that the rhs of Eq. ( |4^ ) is the generating function 

+00 

nni,n2;e) = J]i^i(ni,n2)e^ C < 1, (46) 
i=o 

of Pj (721,712) - the probability that a particle performing some special type of random walk 
on the sites of a two-dimensional square lattice and starting at the origin arrives exactly on 
the j-th step to the site with the lattice vector Tiiei -|- n2e2- This special random walk is 
characterized by A(k) - the structure factor (cf. ref. [^), which defines a biased random 
walk with the following elementary jump probabilities 



2{,/7[^ + A2)' 



in the directions it ei of the lattice, (47) 
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, in the directions it e2 of the lattice. (48) 



Hence, -Fni,n2 can be interpreted as 

/A ^ \ "1/2 

^ni,n, = i-^] P{m,n2;0, (49) 

and thus can be thought of as the generating function of a two-dimensional biased random 
walk. 

5 Asymptotical behavior of the density profiles at large dis- 
tances from the tracer particle. 

Equations (^) to (^) allow us to deduce an asymptotical behavior of the density profiles in 
front of and past the stationary moving TP. To do this, note first that these asymptotic forms 
of /i„.o with n ±00 can be most straightforwardly obtained from the generating function 
of hnfi-. 

+00 

N{z)= Vo^", (50) 



which we proceed now to calculate. Using Eqs.(32) and (|37D, we find that 



z{z - 1) (Aih{ei) + psiAi - A^i)) +{l-z) (A_i/i(e_i) - - ^_i)) 
Niz) = ^ ^ ^ ^ 

V(Z - Zi){z - Z2){Z - Z^){Z - Za) 



"(e.) , '-"-''i'^'^^! -! ■ (51) 



^ (z - zi){z - Z4) 
where the roots 2,- are defined as 



.1 = ^ f + ^.-^(f + A,)'-^,A_, , (52) 



1 « A 
Z2 = -. { A2 

A.\2 



a \ 2 

--^sj -AiA^i 



(53) 



1 « A 



a \2 

--Aaj -AiA^i 



(54) 



(55) 
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which obey 

< zi < Z2 < 1 < Z3 < Z4. (56) 

In consequence, N{z) is a holomorphic function in the annular region of inner radius Z2 and 
the outer radius Z3. As explained in Appendix B, the behavior of when n +00 (resp. 
n — > —00) is controled by Z3 (resp. z^^). 

5.1 Asymptotics of density profiles at large distances in front of the tracer 
particle. 

The asymptotical behavior of the density profiles in front of the TP is supported by the 
behavior of the generating function N(z) in the vincinity of Z3 (see Apendix B for more 
details). We find then that for n — > +00 the density obeys 

Z3(z3-l)(^Ai/lei+/9s(^l-^-l)) +(l-^3)(^-l/le-i-/0s(^l-^-l)) 1 1 

A_1V^(Z3 - Zi){z3 - Z2){Z4 - Z3) ^Zl' 

(57) 

which means that in front of the TP the deviation hnfi decays exponentially with the distance, 

exp ( — n/A-i- I 

/i„,o ~ ^^7^ '- , (58) 

where the decay amplitude obeys 

(Z3 - 1) (^i/i(ei) + p,{A^ - ^_i)) + (I/Z3 - 1) fA_i/i(e_i) - p,{A^ - A^^, 



2 



W,^Zi^V_^J . (59) 



-1/4 

while the characteristic decay length is given by 

A+ = l/ln(z3). (60) 
Note that here A-|- is finite for any values of the system parameters. 

5.2 Asymptotics of density profiles at large distances past the tracer par- 
ticle. 

We first note that one of the roots of the generating function, namely Z2, gets equal to 
unity when both / and g are strictly equal to zero, which results in the exact cancellation 



of the multiplier ^Jl — z both in the nominator and the denominator in the first term on 



the rhs of Eq.(51). This signifies that in the limit when exchanges with the reservoir are 
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forbidden, the singular behavior of the generating function at the vicinity of Z2 is essentially 
different compared to the case when both / and g are greater than zero. Consequently, one 
has to consider separately the behavior in the case of non-conserved particles number, when 
exchanges with the reservoir persist, and the case when both / and g are equal to zero while 
their ratio is kept fixed. 



A. Non-conserved particles number. In case when both / and g have non-zero values, 
we find (cf. Appendix B) that when n +00 the density deviation follows 

exp ( — n/\\ 

h.n,o-K^ ^— ^, (61) 

where the amplitude is given by 

{Z2 - 1) (Ai/lei + Ps{Al - A_i)) + (1/Z2 - 1) (^-l/ie_i " Ps{Al - A^i) 




-1/4 

M ) ' <«^' 

and the characteristic decay length obeys 

A_ = -l/ln(z2). (63) 

Consequently, in the case when particles exchanges with the reservoir are permitted, the 
density deviation from the equilibrium value ps decays exponentially with the distance from 
the TP. Note that the decay lenghts satisfy the inequality A_ > A-|-, which means that the 
correlations between the tracer and particles of the monolayer are always stronger past than 
in front of the TP. Note also that A_ diverges when both / and g tend to zero, which signals, 
as we have already remarked, that the decay of hnfl may proceed differently in this case, 
compared to the exponential dependence in Eq. (|6TD . 



B. Conserved particles number. Suppose now that both / and g are equal to zero, 
while their ratio is fixed and obeys f /g = Ps/{^ — Ps)- As we have already remarked, this 
situation corresponds to the customary model of a two-dimensional hard-core lattice gas 
without exchanges with a reservoir. For this situation, we find (see Appendix B for more 
details), that when n +00 the deviation of the particle density from the equilibrium value 
Ps obeys 

= (1 + 1; + fQ), (64) 
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where the decay amphtude K'_ is given by 



K'_ = -_<^ (Ai/i(ei)-A_i/i(e_i) + 2p,(^i-^_i)) 



1 



+Me.)^^^} (65) 

Remarkably enough, in this case the correlations between the TP position and the particles 
distribution vanish algebraically slow with the distance! This implies, in turn, that in the 
conserved particles number case, the mixing of the monolayer is not efficient enough to prevent 
the appearence of the quasi-long-range order and the medium "remembers" the passage of the 
TP on a long time and space scale, which signifies very strong memory effects. We note also 



that the algebraic decay of correlations in this model has been predicted earlier in Ref.|20|. 
However, the decay exponent has been erroneously suggested to be equal to 1/2, as opposed 
to the value 3/2 given by Eq.(|6^). As well, the amplitude K'_ happens to have a different 



sign, compared to that obtained in Ref.|20|, which invalidates the conclusion that the overall 
relaxation to the equilibrium value ps might show a non-monotoneous behavior as a function 
of the distance past the TP. 



6 Formal expression for the density profiles and general force- 
velocity relation. 

The Eqs.(^9|) and ([42|) display /ini,n2 &s a function of the coefficients Ay that remain to be 
determined. As a matter of fact, these coefficients depend themselves on the local densities 
in the immediate vicintiy of the tracer, i.e. on h{ey). This implies that we have to determine 
them from Eqs.(p9|) and (^) in a self-consistent way. 

Setting in Eq.(H) A = ejy, where u = {±1, 2}, results in the following system of equations 



where 



and 



C 



Ch 


= Ps{Ai 


-A. 


( /i(ei) \ 






/i(e-i) 


, F = 


(Vi 


V Ke^) ) 







/AiV^iF,, - 

^lV_iFe_, 
V ^lV_iFe, 



-V-l)Fe, 
V-l)Fe_, 
- V_i)Fe, 



a .4_iViFe, 



A2V-2F^_, 



(66) 



(67) 



(68) 
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The solutions of this 3x3 hnear system are 

yu = {±l,2}, h{e^) = ps{A^-A_{)^^, (69) 

det C 

where Cy stands for the matrix obtained from C by replacing the i/-th column by the column- 
vector F. The substitution of expression in Eq.(p9D into the definition of the coefficients A^,, 
([T^ ) , results in the following system of three equations 

Vz. = {±1, 2}, A, = l + —p, ll-ps- PsiA, - ^-i)^^ I , (70) 

T [ detC J 

which implicitly determines all unknown coefficients Ai, and hence, the local deviations /i(e^), 
defined as 

h{e^) = {l-p,) + -^{l-A,). (71) 

Finally, for every lattice point, except for the origin occupied by the TP, the density profiles 
are given by Eq.(3S), where -Fni,n2 is defined by Eq. (^). 

Next, substituting Eq.([69|) to Eqs.(^) and (|T9[), we find that the TP terminal velocity 
obeys 

Vtr = -\ [Pi -P-i)(l - Ps) - Ps{Ai - A_i) \, (72) 

T (■ det C J 

which can be rewritten, taking into account Eq.(p^, as 

T/ u-i ^/^ , 4T*pidetC'i -p-idetC-n -1 

Vtr = -{pi-p-i)[l- ps)<l + Ps —~ \ ■ (73 

T I r detC J 

This last equation represents the desired general force- velocity relation for the system under 
study, which is valid for arbitrary magnitude of external bias and arbitrary values of other 
system's parameters. In the general form as it is, Eq.([73|) is apparently not very useful, but 
may describe physically interesting behavior in two particular limits - the case of infinitely 
strong bias, when the TP performs totally directed random walk, which is appropriate to 
motion of the AFM tip, and the case when the bias is vanishingly small. While the former 



case has been studied in detail in Ref. |Q , the results for the latter case have been only briefly 
outlined in Ref.|^^. In the next Section we present detailed derivation of the limiting form 
of the force-velocity relation in case of a vanishingly small external force. 



7 Vanishingly small external bias. 

We turn now to the limit PE <C 1, in which case the problem simplifies considerably and 
allows to obtain explicit results for the local densities in the immediate vicinity of the TP 
and consequently, the TP terminal velocity and diffusivity. 
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7.1 Local particle density in the vicinity of the tracer in the limit (3E <^ 1. 

Setting A = ei and A = e_i in Eq.(^9|), we obtain equations obeyed by h{e±i): 

i/i(ei) = J^A,/i(e,)V_,Fe, -^-i)(Vi - V_i)Fe,, (74) 



and 

Q/i(e_i) = J];^,/i(e,)V_,Fe_, - PsiAi - A_i)(Vi - V-i)F,_,. (75) 

Consequently, the discontinuity in the density profile is given by: 

5/1= (/i(ei)-/i(e_i)) = a-i^^,/i(e,)V_,(Fe, -Fe_J- 

- p,{A^ - ^-l)(Vi - V_i)(Fe, - Fe_,) (76) 

Our aim is now to compute the leading order contribution to bh in the limit fiE ^ 1. To this 
purpose, we first expand /i(e±i) and Ay in the Taylor series in powers of E retaining only 
linear with the field terms: 

P. = ^(l + ^(E-e.) + 0(£;2)), (77) 
/lei oc E and = -h^^ + 0{E^), (78) 



/ie±, = OiE^), (79) 



and consequently, 



Ay = 1 + -{1- ps) + -((1 - Ps)^{^ ■ e.) - h{e^)) + 0{E^), (80) 



where we have made use of the symmetry relations given by Eq. (|l5|). Using these equations. 



we can expand the terms in the right-hand-side of Eq.(76) as: 



^Ayh{e^)V^y{F,, - F,_,) = Y,Ayh{e^)V.y{F,,iE = 0)-F,_,iE = 0)) + OiE^ 



/ie,(AiV_i - A_iVi)(Fe,(F = 0) - F,_,{E = 0)) + 0{E 
2/i(ei) ^ ^ 2 



2^ 



l + -{l-p,)\+OiE'), (81) 



C{2Ao/ao) V r 

where 

fOO , ^ —1 

10 



C{x) = l e"*((Io(xt)-l2(xt))Io(xt)dtl , (82) 
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Ao = A,{E = Q) = l + —{l-ps), (83) 

T 



and aQ = a{E = {}). Next, we have that 

((1 - ps)(3aE - 2/i(ei)) + 0{E^) (84) 



p,(Ai-A_i)(Vi- V_i)(Fe, - i^e_J = 

2ps 



£(2Ao/ao) r 

Substituting Eqs.(|8l|) and (|8|) into Eq.(^, we obtain 



aoi'(2Ao/ao) V / 
+ .^9^; / . -((l-Ps)/?a£;-2fe(ei)) + 0(i^^), (85) 

which leads to the desired result for the deviation of the local density just in front of the TP 
from the equilibrium value ps- The leading order contribution to h{e\) with respect to the 
field is thus given explicitly by 

h{e^) = Ps{l- Ps)l3aE—[aoC{2Ao/ao)-l + {2,ps-l) — \ +0{E^). (86) 



Hence, the discontinuity 6h in the monolayer particles density in the immediate vicinity of 
the TP equals, by virtue of Eq.([78[), twice the expression on the right-hand-side of Eq. (|8^) . 

7.2 Friction coefficient and the Stokes formula. 

Expanding the TP velocity Vtr in the Taylor series in powers of E and again, retaining only 
linear with E terms, we get 

Vtr = - Ps)E - ^ (/i(ei) - /i(e_i)) + 0{E^) (87) 



Next, making use of Eqs.(|78|) and (SC), we arrive at the following explicit result 



Vtr ~ -;— (1 - Ps)E< 1 



4r \ r ao/:(2Ao/Qo) - 1 + (3p, - 1)t*/t 

which signifies that in the limit of a vanishingly small external bias the frictional force exerted 
on the TP by the monolayer particles is viscous. Note also that Eq.(|88[) is quite similar to the 
well-known Stokes formula, i.e. Vtr ~ E/C,, and in our case the friction coefficient C, is given 
explicitly by 



/3(j2(l -ps)\ T ao/:(2^o/ao) - . 
Note now that the friction coefficient ^ is the sum of two contributions. The first one, 

/3CJ^(1 - p,) 
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is a typical mean-field-type result and corresponds to a perfectly homogeneous monolayer (see 
discussion following Eq.(6)). The second one, 

C = ""-Pi (91) 

has, however, a more complicated origin. Namely, it is associated with the cooperative behav- 
ior emerging in the monolayer - dehomogenization of the particle distribution in the adsorbed 
monolayer due to the presence of a driven impurity (the TP) and formation of stationary 
density profiles, whose characteristic properties depend on the velocity Vtr- 

7.3 Diffusion coefRcient and Einstein relation 

Lastly, assuming a priori that the Einstein relation holds for the system under study, we 
estimate the TP diffusion coefficient Dtr as 

Dtr = r'c' = ^(1 -Ps)\i- , , — j^rr] ■ (92) 

4r [ r aoL{2Ao/ao) - 1 + [Sps - 1)t*/t J 



It seems now interesting to compare our general result in Eg .(19 21) against the classical result 
of Nakazato and Kitahara [l^, which describes TP diffusion coefficient in a two-dimensional 
lattice gas with conserved particles number. Setting / and g equal to zero, while assuming 
that their ratio has a fixed value, f /g = Ps/{^ ~ Ps)i "we have then that 

Noticing next that 

= ^hm (P(0, 0; - P(2, 0; o) , (94) 
where P(ni,n2;^) has been defined by Eq.(E^, and using the fact that [p|] 



hm P(0,0;0-mO;e) =4--, (95) 

we find that the right-hand-side of Eq.(^) attains the form 

- _ J 2psT* l-2/7r \ / \ 

A. - ^ (1 - P.) 1 1 - ^ 1 + (1 _ p^)^*,^ _ (1 _ 2/vr)(l + (1 - 3p.)rVr) / ' ^^^^ 



which expression is exactly the same as the one obtained earlier in Refs.[^ and [17| within 
the framework of different, compared to our approach, analytical techniques. The result in 
Eq.(|9^) is known to be exact in the limits *S 1 and ~ 1, and serves as a very good 
approximation for the self-diffusion coefficient in hard-core lattice gases of arbitrary density 

HI. 
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8 Conclusion 



To conclude, we have studied analytically the intrinsic frictional properties of 2D adsorbed 
monolayers, composed of mobile hard-core particles undergoing continuous exchanged with 
the vapor. Our analytical approach has been based on the master equation, describing the 
time evolution of the system, which has allowed us to evaluate a system of coupled dynamical 
equations for the tracer particle velocity and a hierarchy of correlation functions. To solve 
these coupled equations, we have invoked an approximate closure scheme based on the decom- 
position of the third-order correlation functions into a product of pairwise correlations, which 
has been first introduced in Ref.[Q for a related model of a driven tracer particle dynamics 
in a one-dimensional lattice gas with conserved particles number. Within the framework of 
this approximation, we have derived a system of coupled, discrete-space equations describing 
evolution of the density profiles in the adsorbed monolayer, as seen from the moving tracer, 
and its velocity Vtr- We have shown that the density profile around the tracer is strongly 
inhomogeneous: the local density of the adsorbed particles in front of the tracer is higher 
than the average and approaches the average value as an exponential function of the distance 
from the tracer. On the other hand, past the tracer the local density is always lower than the 
average, and depending on whether the number of particles is explicitly conserved or not, the 
local density past the tracer may tend to the average value either as an exponential or even 
as an algebraic function of the distance. The latter reveals especially strong memory effects 
and strong correlations between the particle distribution in the environment and the carrier 
position. Next, we have derived a general force-velocity relation, which defines the terminal 
velocity of the tracer particle for arbitrary applied fields and arbitrary values of other system 
parameters. We have demonstrated next that in the limit of a vanishingly small external bias 
this relation attains a simple, but physically meaningful form of the Stokes formula, which 
signifies that in this limit the frictional force exerted on the tracer by the adsorbed monolayer 
particles is viscous. Corresponding friction coefficient has been also explicitly determined. In 
addition, we estimated the self-diffusion coefficient of the tracer in the absence of the field 
and showed that it reduces to the well-know result of Refs.||l^ and jl^ in the limit when the 
particles number is conserved. 
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Figure Captions. 

Fig.l. Adsorbed monolayer in contact with a vapor. Grey spheres denote the monolayer 
(vapor) particles; the grey sphere with an arrow stands for the driven tracer particle. 

Fig. 2. Singular points Zi of the generating function N{z) and the integration contour 
(thick line). 

Fig.3. Piecewise contour H{n) of the Hankel type. 
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Appendix A 



In this Appendix we present the details of calculation of different contributions to the 
time evolution of the pair correlation function k{\;t) (cf. Eq.(|^), associated with four lines 
in the master equation @ . 

Consider first the contribution associated with the hopping motion of the adsorbed par- 
ticles: 

CiW ^ E E E r7(Rtr + A){P(Rtr,r/'-'^;t)-P(Rtr,r/;t)} 

= E E E {??'"''^(Rtr + A)-7?(Rtr + A)}P(Rtr,r?;t). 

R-tr/'; M=l)2 r^^^Rtr-e^.Rtr 

(97) 

Now, one has then to distinguish between two possible situations: 

(a) when A 7^ ej^, one has 

^ {r?'-''^(Rtr + A)-r/(Rtr + A)}P(Rtr,?7;i) 

^J.=1,2 r^Rtr-e^i.Rtr 
= E + ^-M^tr + A) 

= J^V^r?(Rtr + A), (98) 
which yields, eventually, 

Ci(A) = J]V^A;(A;i). (99) 

(b) when A = ej^, i.e. at the sites adjacent to the tracer, we find 

E E {??"'^(Rtr + e^)-r/(Rtr + e^)}P(Rtr,r/;t) 

/^=l,2 r^Rtr-e^,Rtr 

= (E^'^-^-'')'?(^tr + e^), (100) 
which implies that 

(^ilet.) = (E^M-^--)^(e«^;*)- (101) 

Finally, using the Kroneker-delta 5^,6^ , we can generalize both results for arbitrary A: 

Ci(A) = (E - Vm^-m)^!^; t)- (102) 
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Next, we turn to the contibution stemming out of random biased hopping motion of the 
tracer particle. This reads 

C2(A) = — ^ ^p^r/(Rt, + e^){(l-7?(Rtr))P(Rtr-e^,ry;f) 

- {1 - r](Rtr + e^)) P(Rtr,v;t)} 
= ^ - r?(Rtr + e^)){??(Rtr. + + e^) - r/(Rt,. + ei,)}P(Rtr,r/;t) 

Rtr,»? M 

= ^ 5^PM(l-?7(R.tr + e^))V^?7(Rtr + A)P(Rtr,r/;t), (103) 

Rtr,»? M 

Further on, we consider the contribution associated with desorption of the monolayer 
particles: 

CsiX) ^ ^gY^Yl E^(^tr + A){(l-r?(r))P(Rtr,77'';i)-^(r)P(Rtr,r?;t)} 
= ^^^E E E^Wi^'(^tr + A)-??(Rtr + A)}P(Rtr,r?;t) 

Rtr r^Rtr V 

= 49^^7?(Rtr + A)(l-277(Rtr + A))P(Rtr,?7;t) (104) 

Rtr V 

Taking into account that ?7(R) assumes only two values - and 1, and hence, that r/^(R) = 
rj(R}, we have 

CsiX) = -igkiX-t). (105) 

Lastly, the contribution due to adsorption of the particles from the vapor phase onto the 
lattice reads 

C4(A) = 4/^ 7?(Rtr + A){7?(r)P(Rtr,r?'-;t)-(l-r?(r))P(Rtr,r?;t)} 

Rtr,r) FT^Rtr 

= ^fY. E (l-^(r)){??'"(Rtr + A)-ry(Rtr + A)}P(Rtr,?7;0 

Rtr,r) FT^Rtr 

= 4/ ^ (l-r/(Rtr + A))(l-27?(Rtr + A))P(Rtr,r/;t) (106) 

Rtr,r) 

Again, since ^/^(R) = ^/(R), the latter equation reduces to 

C4(A) = -4/fc(A;t) + 4/. (107) 
Summing up all four contributions, we arrive at the evolution equation in Eq.(^). 
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Appendix B 



We present here a derivation of the asymptotic behavior of the density profiles at large 
separations from the stationary moving tracer. The generating function N{z), defined by 
Eq.dH), has the roots Zi satisfy Eqs.(|52|) to (|55|). Introducing the cuts depicted on Fig. 2, 



one observes that N{z) is an analytic function in the annular region of inner radius Z2 and 
the outer radius 2:3. In consequence, density deviations hnfl from the equilibrium value ps at 
positions A = ein are given by the Cauchy formula 

'-°=2^/^^^^' ^'''^ 

where the contour of the integration is a positively oriented circle, centered around O, of 
radius R such that Z2 < R < Z3 (see Fig. 2). 

I. Asymptotic behavior of h^fl in the limit n — > +00. Following closely the reasonings 
of Flajolet et al. |3^, we notice that the asymptotical behavior of hnfi in the limit n — > +00 
is supported by the behavior of the generating function N{z) in the vicinity oi z = Zi, where 



Zi is the first root encountered when one tries to deforme the contour in the integral (108) by 
increasing its radius. In our case, the relevant root is z^. Now, we have to choose a contour 
of integration that comes close enough to "capture" the behavior of N{z) in the vicinity of 
this leading singularity. Let us formalize this idea. 

We begin by expanding N{z) in the vicinity of Z3, which gives 

N{z) = C{zs - - h{e2) + 0(^/i^), (109) 

where C is a constant defined by 

^ Z3{z3 - 1) (^i/i(ei) + ps{Ai - A^i) ) + (1 - Z3) ( ^_i/i(e„i) - p,{Ai - A^i) ' 



'a/2-A2Y Al 



V8A2Z3 

-1/4 



Then, Eq.(|108|) attains the form 



(110) 



,„„ = ^ife^,,_MSH)/ 1 ,.+i/5(4Ild.. (Ill) 



=0 

Consider next the first term on the rhs of Eq.(|lllD, i.e.. 
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One notices first that the contour of integration in Eq.(I12) may be replaced by the piecewise 
contour Ti.{n) of the Hankel type (cf. Fig.3): 



where 



n{n) = n-{n) + n+{n) + n^{n), 



n-{n) = {z = w- ^, w>l} 
n+{n) ={z = w+^, w>l} 



(113) 



(114) 



Changing next the variable of integration in Eq.(|lll|), z = 23(1 + t/n), we have that Dn obeys 



Further on, expanding the kernel 

, \ — n— 1 
1 + - 

n 



1/2 -n^^ 



-1/2 



2m 



n 



i-ty 



-1/2 



1 + 



n 



-n— 1 



dt 



-(n+l)ln(l+t/") = e"*(^l + 



f-2t 



2n 



+ 



(115) 



(116) 



we observe that the integrand in Eq. (|115| ) converges to (— t)~^/^e~*, which is just the kernel 
appearing in the Hankel's representation of the Gamma function. Turning to the limit n — > 
+00, and using an approximation 



1 + 



t 



n 



-n-l 



1 + 



(117) 



in the integral in Eq.( 115 ), we find eventually that 

C 1 



^r(i/2) 

C 1 



l + O 



n 



l + O 



n 



yfz^ ^z^ 

Next, using essentially the same kind of arguments, it is easy now to show that 



1 /0(V?^) 



2i7r 



dz = 



1 



3/2 



n 



which implies that in the limit n — > +00, the leading behavior is given by 

C 1 



n,0 



(118) 



(119) 



(120) 



which is exactly the result in Eqs.(pq) to (pO 
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II. Asymptotical behavior of /i„ o in the Umit n — > — cxd. To study asymptotical be- 
havior of tiie density profiles at large separations past stationary moving tracer, we proceed 
exactly in the same way as in the previous paragraph. First, we introduce the generating 
function of the form 

+00 

M{z)= h-nfiz'' (121) 

n=— 00 

Since M{z) obeys M{z) = N[l/z), one immediately has that M{z) is given explicitly by 

(1 - z) (^i/i(ei) + p,{Ai - A^i)) + z{z - 1) f^_i/i(e_i) - ps{Ai - A^i 

M{z) = ^ ^ ^ 

Ai^{z - z^^){z - z:,^){z - z^^){z - Z4-I) 

+ K..)iJf^^^^-l]. (122) 
\\l {z - z^ )[z - z^ ) ) 

Hence, M[z) is an analytic function in the annular region centered around O, of inner radius 
z^^ and the outer radius Z2^ . 

Now, the leading singularity for M{z) is in the vicinity of 2: = z^^. As explained in 
the text, the nature of this singularity Z2^ is different depending on whether the number of 
particles in the monolayer is explicitly conserved or not. 

A. Non conserved particles number. In this case Zg^^ 7^ 1 and hence, in the vicinity 
of z = Z2^ ^ the generating function M{z) behaves as 



M{z) = E{z2^ - z)-i/2 _ ^(g^) ^ 0{^Z2^ - z), (123) 
where £^ is a constant, defined by 

(1 - 2,-1) (Ai/i(ei) + p.iAi - A^i)) + Z2\z2^ - 1) (^_i/i(e_i) - ps{Ai - A_i) 
E = ^ 



a/2-^2^^ A-i 



A, / A, 



-1/4 

(124) 



Inverting Eq.( |123| ), we find then the following asymptotical result: 



EJ^ 11 , ^ 

h-n,o-^^-^, (125) 
\/\n\ ^2 



which corresponds to our Eqs.(|6l|) to (|63[). 
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B. Conserved particles number. In this case, z^^ = ^ and consequently, one has 
that M(z) admits the fohowing expansion in the vicinity oi z = z^^, 

M{z) = -/i(e2) + G(l - + O ((1 - z)3/2) , (126) 

where G is the constant given by 



1 [Aih{er)+ps{Ai-A_i))-[A^ih{e^r)-ps{Ai-A^i] 

(127) 

Inverting Eq.( |12£| ), we find that in the conserved particles number case, the asymptotical 
behavior of the density profile at large separation past the stationary moving tracer particle 
obeys 

i.e. the behavior described by our Eqs.(|64[) and (pBI). 
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Figure 1: PRB, Benichou et al. 




Figure 2: PRB, Benichou et al 
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Figure 3: PRB, Benichou et al 
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